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Summary. — Dusty hyperluminous galaxies in the early universe provide unique 
environments for studying the role of massive stars in the formation and destruction 
of dust. At redshifts above 6, when the universe was less than 1 Gyr old, dust could 
have only condensed in the explosive ejecta of Type II supernovae (SNe), since most 
of the progenitors of the AGB stars, the major alternative source of interstellar 
dust, did not have time to evolve off the main sequence. We present analytical 
models for the evolution of the gas, dust, and metals in high redshift galaxies, with 
a special application to SDSS J1148+5251, a hyperluminous quasar ai z — 6.4 
hereafter referred to as J1148+5251. We show that an average SN must condense 
at least 1 Mq of dust to account for the mass of dust in this object, when grain 
destruction by supernova remnants (SNRs) is taken into account. This required 
yield is in excess of ^ 0.05 Mq, the largest mass of dust inferred from infrared 
observations of Gas A. If the yield of Gas A is typical, then other processes, such 
as accretion onto preexisting grains in molecular clouds are needed to produce the 
mass of dust in J1148+5251. For such process to be effective, SNR must significantly 
increase, presumably by non-evaporative grain-grain collisions during the late stages 
of their evolution, the number of nucleation centers onto which refractory elements 
can condense in molecular clouds. 
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1. The presence of massive amounts of dust at high redshift 

The detection of massive amounts of dust in hyperluminous infrared (IR) galaxies at 
redshifts z > 6 raises challenging questions about the sources capable of producing such 
large amount of dust during the relatively short lifetime of these galaxies (e.g. [1, 2]). 
For example, the galaxy SDSS J1148-)-5251 (hereafter J1148-F5251) located at z = 6.4 
was observed at far-IR and submillimeter wavelength [3, 4, 2]. Figure 1 depicts the 
observed fluxes and spectral fits of various dust emission models through the data. The 
average IR luminosity of the source is Lj^ ^ 2 x 10^^ Lq, and the average dust mass 
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is Md ~ 2 X 10* Mq. Using the Kennicutt relation [5], one can derive a star formation 
rate (SFR) of ~ 3000 Mq yr~^ from the observed far-IR luminosity. For comparison, 

the Milky Way galaxy is about 10 Gyr old, has an average SFR of ^ 3 Mq yr~^, and 
contains about 5 x 10'' Mq of dust, a significant fraction of which was produced in AGB 
stars. 
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Fig. 1. - Spectral fits of several possible dust compositions to the observed far-IR and submil- 
limeter observations of J1148+5251. Wavelengths are in the observer's frame of reference. Dust 

masses range from (1 — 5) x 10** A4q and luminosities range from (2 — 3) x 10^^ Lq, depending 
on dust composition. References to the observations can be found in [6]. 



At z = 6.4 the universe was only 890 Myr old, using standard ACDM parameters 
{n,n = 0.27, Qa ^ 0.73, and Hq = 70 km s'^ Mpc"!). If Jl 148+5251 formed at z = 10 
then the galaxy is only 400 Myr old. If the SFR had occured at a constant rate over 
the lifetime of the galaxy, its initial mass should have been about 10^^ Mq, which is 
significantly larger than the dynamical mass Mdyn « 5 X 10^° Mq of the galaxy [7]. The 
high observed SFR may therefore represent a recent burst of star formation that has 
lasted for only about 20 Myr. The galaxy J1148+5251 is therefore at most ~ 400 Myr 
old, and probably significantly younger with an age of only ^ 20 Myr. Adopting a 
current gas mass of Mg = 3 x 10^° Mq for this galaxy we get that the gas mass fraction 
at 400 Myr is about 0.60. The dust-to-gas mass ratio is given by Zd = Md/Mg = 0.0067. 

A significant fraction of the dust in the Milky Way was produced in AGB stars. 
However, these stars are not likely to contribute significantly to the formation of dust 
in very young galaxies, since the low mass stars (M w 2 Mq) that produce most of 
the dust did not have time to evolve off the main sequence [8, 9, 6]. In contrast, core 
collapse supernovae (SNe) {M > 8 Mq) inject their nucleosynthetic products back into 
the ISM shortly {t < 20 Myr) after their formation, resulting in the rapid enrichment of 
the interstellar medium (ISM) with the dust that formed in their explosive ejecta. But 
can SNe account for the large amount of dust seen in this object? The answer to this 
question is complicated by the fact that SNe are also the main source of grain destruction 
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during the remnant phase of their evolution [10, 11]. The problem can therefore only be 
addressed with a detailed model for the evolution of dust in these systems. 

2. — The evolution of dust in evolved galaxies 

The evolution of dust is best studied in the framework of chemical evolution models 
that follow the galactic metal enrichment resulting from stellar nucleosynthesis (e.g. [12]). 
Such chemical evolution models need to be generalized to include processes unique to the 
evolution of dust: the condensation efficiency of the refractory elements in stellar ejecta, 
the destruction of grains in the ISM by expanding supernova remnants (SNRs), and the 
growth and coagulation of grains in clouds [8]. 




Fig. 2. - The evolution of SN-condensed silicate and carbon dust and AGB-condenscd carbon 
dust as a function of time. The model used to produce the figure is described in the text, and 
is presented in more detail in [13]. 



Figure 2 shows the evolution of SN and AGB-condenscd dust as a function of time. 
Details of the dust evolution model are given in [13]. SNe can form both silicate and 
carbon dust since their ejecta have layers with different C/0 ratios that may be larger 
or smaller than unity which are never microscopically-mixed. In contrast, the ejecta of 
AGB stars are either carbon- or oxygen-rich. Ideally, AGE stars are therefore expected 
to form either carbon or silicate dust. Observations show that some stellar outflows 
may have more complex chemistry which may be the results of mass loss from a binary 
system, since they show the presence of both silicates and C-rich dust products [14]. For 
simplicity, the evolutionary models used in Figure 2 assume that stars with masses below 
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8 Mq form either carbon or silicate dust depending on the C/0 ratio in their ejecta, and 
that higher mass stars that end as core-collapse SNe form both silicate and carbon dust 
according to the prescription used by [8]. For AGB stars, the figure depicts only their 
contribution to the carbon dust enrichment of the ISM. The most massive AGB stars 
are 0-rich, so that the onset of the AGB contribution to the carbon dust abundance is 
delayed until the C-rich stars evolve off the main sequence. However, the IMF-weighted 
yield of carbon dust peaks at stellar masses of about 2 Mq, so most of the carbon dust 
is injected into the ISM with a delay time between 0.4 and 4 Gyr, depending on the star 
formation history and the grain destruction rate. The subsequent decline in the yield of 
carbon dust results from the fact that lower mass AGB stars are 0-rich and will only 
form silicates in their ejecta. 




Fig. 3. - Comparison between the metallicity trends of the PAH abundance derived from the 
observed SED and those derived from the chemical evolution model. The shaded area represent 
the range of model prediction for different grain destruction and star formation rates. Details 
of the figure are described in [13]. 
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The qualitative features of this dust evolution model were recently confirmed by 
mid-IR photometric and spectral observations of the features of polycyclic aromatic hy- 
drocarbon (PAH) molecules in nearby galaxies with the Spitzer and ISO Infrared Space 
Observatory satellites (see [13] for references). The observations showed a deficiency in 
PAH emission from low metallicity galaxies, and a trend of increased emission in galax- 
ies with higher metallicity. The PAH abundance, derived from detailed modeling of the 
spectral energy distribution of nearby galaxies [13], also correlated with their metallic- 
ity. This trend can be explained as a natural consequence of stellar evolutionary effects. 
Since PAHs are made in AGB stars, the trend simply reflects the delayed injection of 
these macromolecules into the ISM [15, 13]. Figure 3 is a mapping of the trend of PAH 
abundance with age onto a trend with metallicity for the same chemical evolution model. 
The figure illustrates the disctinct evolutionary trends of SN- and AGB-condensed dust. 
The shaded area represents the range of model predictions for different grain destruction 
and star formation rates. Details of the chemical evolution model can be found in [13]. 



3. — The chemical evolution of dust in young systems 

The detailed dust evolution models described above show that the contribution of 
AGB stars to metal and dust abundance can be neglected in galaxies with ages less than 
about 400 Myr. The equations for the chemical evolution of the galaxy can then be 
solved using the instantaneous recycling approximation, which assumes that stars return 
their ejecta back to the ISM promply after their formation. In the following wc present 
the results for a closed box model. The case for the infall model is given in [6] . 

In the closed box model, the mass of the ISM changes due to the process of star 
formation (astration) at a rate given by: 

(1) ^ = _(i_ij)^(t) 

where R is the fraction of the stellar mass that is returned back to the ISM during the 
stellar lifetime. Assuming that the SFR, 'ip{t), is proportional to the gas mass to some 
power k, we can write that: 



where Mq is the mass of the system which initially consists of only gas, Hg{t) is the 
fractional mass of gas at later times t, and tpo is a parameter that is adjusted to fit the 
observed gas mass fraction a,t t = 400 Myr, taken to be the age of the galaxy. 

The equation for dMg/dt can then be written in terms of the gas mass fraction as: 



(3) 



Analytical solutions are available for arbitrary values of k: 



(4) 



IJ.g{t) = exp 



k = l 
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with the initial condition that Hg = 1 at time t = 0. 

The equation for the evolution of the mass of dust, in the galaxy is given by: 



(5) 



dt 



-Zd ilj{t)+YaRsN{t) 



Md{t) 



The first term represents the destruction rate of dust by astration, the second term its 
formation rate by SNc, and the third term the destruction rate of the dust by expanding 
SN blast waves. In the equation above, = M^/Mg is the dust-to-gas mass ratio, 
is the average yield of dust in Type II supernovae, Rsn is the galaxy's supernova rate 
determined by the SFR and the stellar IMF, and Td is the lifetime of the dust against 
destruction by SNRs. The lifetime for grain destruction is given by: 



(6) 



Td 



Mg{t) 



(misM) Rs 



N 



where (jtiism) is the effective mass of ISM gas in which the dust grains are completely 
destroyed by a single SNR. In our Galaxy, Mg « 5 x 10^ Mq, Rsn ~ 0.02 yr-^ and 
{misM) ~ 500 Mq, giving an average dust lifetime of w 5 x 10^ yr [10]. The SN rate can 
be written as: 



(7) 



-Rsn = 



where is the mass of all stars born per supernova event. For example, for a Salpeter 
IMF in the 0.1 and 100 Mq range of stellar masses, m-^— 147 Mq. 

Using equations (3) and (5), the evolution of the dust mass can be expressed as a 
function of the remaining gas fraction as: 



(8) 



MdilJig) = Yd 



_ (toism) +Rrn^ 
where is a dimensionless parameter: 



(9) 



{1-R)m^ 



The expression for Md{i^g) is independent of the star formation history of the galax;y, 
but depends on the details of the IMF. 

Equation (8) can be rewritten to give the dust yield required to obtain a given dust- 
to-gas mass ratio, Zd = Md/Mg, when the galaxy reaches a given gas mass fraction 
lig-. 



(10) 



Yd = Zd 



(wism) + Rfrii, 



1 



Figure 4 shows how much dust an average SN must produce in order to give rise 
to a given dust-to-gas mass ratio, for various grain destruction efficiencies. The value 
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of Yd was calculated when jjg reaches a value of 0.60, the adopted gas mass fraction 
of Jl 148+5251 at 400 Myr. Calculations were performed for two different functional 
forms of the stellar IMF: a Salpctcr IMF in which (f>{in) ^ m""^-^^ and 0.1 < m{MQ) 
< 100; and a top heavy IMF characterized by the same mass limits but a flatter slope 
(f){m) ~ m~^'^°. Here, (p{m) is the number of stars per unit mass interval, normalized to 
unity between 0.1 and 100 Mq. 

The figure shows that, for example, to produce a value of = 0.0067 at /Ug = 0.60, 
a SN must produce about 0.4 (1.2) Mq of dust for a top-heavy (Salpeter) IMF, provided 
the dust is not destroyed in the ISM, that is, (toism) = 0. Even with modest amount of 
grain destruction, (toism) = 100 Mq, the required SN dust yield is dramatically increased 
to about 1 — 2 Mq, depending on the IMF. The horizontal line in the figure corresponds 
to a value of Yd = 0.054 Mq, the largest mass of dust inferred to be present in the ejecta 
of the supernova remnant Cas A [16]. Contrary to the claim by [16], this yield is not 
sufficient to account for the large amount of dust observed in high redshift galaxies, since 
the chemical evolution models of [9] do not include the effect of grain destruction. The 
figure shows that even without grain destruction, the largest observed yield can only give 
rise to a dust-to-gas mass ratio of ~ 4 x 10~^. If the mass of dust in the ejecta of Cas A 
represents a typical SN yield, then other processes, such as accretion onto preexisting 
grains in molecular clouds is needed to produce the mass of dust in J1148+5251. 

4. — Summary 

The early universe is a unique environment for studying the role of massive stars in 
the formation and destruction of dust. The equations describing their chemical evolution 
can be greatly simplified by using the instantaneous recycling approximation, and by 
neglecting the delayed contribution of low mass stars to the metal and dust abundance 
of the ISM. Neglecting any accretion of metals onto pre-existing dust in the interstellar 
medium, the evolution of the dust is then primarily determined by the condensation 
efficiency of refractory elements in the ejecta of Type II supernovae, and the destruction 
efficiency of dust by SN blast waves. 

We applied our general results to J1148+5251, a dusty, hyperluminous quasar at 
redshift z — 6.4 and found that the formation of a dust mass fraction of Zd = 0.0067 in a 
galaxy with an ISM mass of 3 x 10^° Mq, requires an average SN to produce between 0.5 
and 1 Mq of dust if there was no grain destruction (Fig. 4). Such large amount of dust 
can be produced if if the condensation efficiency in SNe is about unity. Observationally, 
the required dust yield is in excess of the largest amount of dust (~ 0.054 Mq) observed 
so far to have formed in a SN. This suggests that accretion in the ISM may play an 
important role in the growth of dust mass. For this process to be effective, SNR must 
significantly increase, presumably by non-evaporative grain-grain collisions during the 
late stages of their evolution, the number of nucleation centers onto which refractory 
elements can condense in molecular clouds. 
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Fig. 4. - The IMF-averaged yield of dust by type II supernova, Yd, that is required to account 
for a given dust-to-gas meiss ratio Zd, is presented for different values of (misM) given in units 
of Mq. Solid and dashed lines correspond to calculations done for a top-heavy and a Salpeter 
IMF, respectively. The horizontal dashed line near the bottom of the figure corresponds a value 
of Yd = 0.054 A/0, the highest inferred yield of dust in supernova ojccta to date [16]. The 
vertical dotted line represents the value of Zd at /Hg = 0.60. Curves are labeled by (misM) given 
in units of Mq . The top two dashed (solid) horizontal lines represent IMF-averaged theoretical 
dust yields for a Salpeter (top-heavy) IMF. 
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